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ABSTRACT 

When  a  radar  guided  missile  engages  a  target  at  low  altitudes  over  the  sea,  the 
radar  signal  reflected  from  the  target  to  the  missile  may  be  subject  to  fading 
owing  to  multipath  effects.  This  fading  can  cause  serious  missile  guidance 
problems.  A  radar  multipath  model  has  been  developed  in  MATRIXx  to 
evaluate  the  performance  of  monopulse  semi-active  missile  systems  against  low 
altitude  targets.  The  multipath  model  takes  into  account  forward  scattering  over 
a  smooth  sea  and  includes  representation  of  the  monopulse  sum  and  difference 
angle  processing  performed  by  the  missile  receiver.  Simulation  results,  showing 
the  adverse  effects  of  multipath  on  the  seeker  boresight  errors,  are  presented  for 
typical  missile-to-target  engagements. 
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THE  PERFORMANCE  OF  SEMI-ACTIVE  RADAR-GUIDED 
MISSILES  AGAINST  SEA  SKIMMING  TARGETS 

EXECUTIVE  SUMMARY 


The  engagement  of  a  sea  skimming  target  by  a  radar-guided  missile  may  be 
subject  to  multipath  effects.  That  is,  the  radar  signal  reflected  from  the  target 
may  return  to  the  missile  seeker  directly  or  via  reflection  from  the  sea  surface 
giving  rise  to  interference  effects.  Such  interference  may  cause  the  target  return 
signal  to  fade.  This  multipath-induced  fading  can  cause  serious  missile  guidance 
problems. 

In  this  report,  a  radar  multipath  return  model  is  developed  in  MATRIXx.  This 
model  is  used  to  evaluate  the  performance  of  a  monopulse  semi-active  missile 
system  against  a  sea  skimming  target.  The  multipath  model  takes  into  account 
forward  scattering  over  a  smooth  sea  and  includes  a  representation  of  the 
monopulse  sum  and  difference  angle  processing  performed  by  the  missile 
receiver. 

A  number  of  parameters  affect  the  return  signal  from  the  target.  These  include 
the  height  of  the  target  above  the  sea  surface,  the  radar  cross-section  of  the 
target,  and  the  sea  state.  Simulations  were  performed  for  a  range  of  these 
parameters  in  order  to  examine  their  impact  on  the  target  returns  and 
subsequently  the  missile  performance. 

The  simulation  results,  which  show  the  adverse  effects  of  multipath  on  the 
seeker  boresight  errors,  are  presented  for  typical  missile-to-target  engagements 
as  a  function  of  the  parameters  identified  above. 

The  multipath  model  described  here  may  be  used  to  assess  the  performance  of 
particular  semi-active  radar  homing  missiles  when  operating  over  a  sea  surface. 
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NOTATION 


Ad  Amplitude  of  the  direct  return  from  the  target 

Ar  Amplitude  of  the  reflected  return  from  the  target 

At  Amplitude  of  the  signal  received  by  the  missile 

G(€)  Antenna  gain  as  a  function  of  the  oflT-boresight  angle  6 

h^  Height  of  the  missile 

hf  Height  of  the  target 

j  Imaginary  number, 

Ji(a:)  Bessel  function  of  the  first  kind  of  order  1 

k  Navigation  Ratio  or  PN  gain 

K  Slope  of  the  monopulse  error  curve 

N'  Effective  guidance  gain 

Nc  Called  for  acceleration 

R  Range  from  the  target  to  the  missile 

Ri  Range  from  the  target  to  the  reflecting  surface 

R2  Range  from  the  reflecting  surface  to  the  missile 

t  Time 

Ta  Autopilot  time  constant 

Tj,  Seeker  gyro  time  constant 

Vc  Closing  velocity 

X  Variable 

A  Difference  voltage  gain 

c  Sea  dielectric  constant 

Complex  dielectric  constant 
a  Boresight  angle  to  the  target  image 

Cm  Boresight  angle  to  the  apparent  target 

€t  Boresight  angle  to  the  target 

7  Turning  rate  of  the  missile  flight  path 

r  Smooth  sea  reflection  coefficient 

A  Radar  wavelength 

Bi  Inertial  angle  to  the  target  image  from  the  missile 

6g  Pointing  angle  of  the  seeker  head  in  inertial  axes 

Bt  Inertial  angle  to  the  target  from  the  missile 

<j>  Sum  of  the  phase  angles  of  the  direct  and  image  returns 

<f>d  Phase  of  the  direct  return  from  the  target 

(f>r  Phase  of  the  image  return  from  the  target 

(l)t  Phase  of  the  signal  received  by  the  missile 

A<^(V’)  Phase  change  on  reflection  for  incidence  angle  ijp 

p  Ratio  of  the  amplitude  of  the  image  signal  to  the  target  signal  at  the  missile 

pc  Coherent  scattering  co-efficient 

Pas  p  modified  to  account  for  the  effect  of  sea  state 


Pi  Ratio  of  image  to  target  amplitudes  before  reflection  of  the  image  signal 

p2  Magnitude  of  the  Fresnel  reflection  co-efRcient 

or  Sea  conductivity 

(Ta  RCS  of  target  as  viewed  by  the  missile  from  a  positive  aspect  angle 

(Tb  RCS  of  target  as  viewed  by  the  missile  from  a  negative  aspect  angle 

&  Turning  rate  of  missile-to-target  line-of-sight  rate 

a  Reconstructed  missile-to-target  line-of-sight  rate 

S  Sum  voltage  gain 

uj  Angular  frequency  of  the  illuminating  radar 

^  Gracing  or  incidence  angle 
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1.  INTRODUCTION 

The  purpose  of  this  report  is  to  examine  the  performance  of  a  generic  radar-guided  missile,  with 
a  monopulse  receiver,  against  low-altitude  targets  flying  over  the  sea  surface.  It  has  been  shown 
that  signal  fading  due  to  the  multipath  interference  can  cause  serious  guidance  problems  for  the 
monopulse  tracker  [1,2].  In  respect  to  a  semi-active  missile,  these  guidance  problems  may  arise  from 
the  multipath  returns  along  any  one  of  the  three  radar  pathways  possible  in  this  type  of  scenario. 
These  pathways  are:  1)  the  transmit  path  from  the  illuminator  to  the  missile  rear  receiver,  2)  the 
transmit  path  from  the  illuminator  to  the  target,  and  3)  the  receive  path  from  the  target  to  the 
missile  front  receiver.  Within  each  of  these  three  distinct  pathways,  there  are  two  possible  forward 
paths  that  the  transmitted  signals  may  take,  namely,  a  direct  path  and  via  reflection  from  the  sea 
surface.  The  multipath  model  presented  here  is  based  on  consideration  of  the  latter  pathway  only 
(see  Figure  1).  Other  multipath-aflfected  pathways  likely  to  pose  serious  guidance  problems  for  a 
semi-active  seeker  will  be  considered  in  future  studies. 

The  modelling  approach  taken  was  twofold.  Firstly,  a  multipath  model  which  adequately  represents 
the  reflection  characteristics  of  the  sea  surface  was  developed  and  implemented  in  MATRIXx.  The 
model  assumes  a  flat  earth  geometry  with  smooth  sea  surface  conditions,  that  is,  zero  sea  state. 
Extensions  to  higher  sea  states  were  modelled  according  to  Ament’s  theory  on  specular  reflection 
[4].  However,  no  modelling  of  the  diffuse  components  of  the  sea  surface  reflections  was  considered 
[3].  The  monopulse  seeker  processing  was  included  in  the  model  as  well  as  attenuation  effects  due 
to  the  sum  and  difference  voltage  gain  patterns  associated  with  the  seeker  antenna.  This  provides 
valuable  insight  into  the  way  the  seeker  processes  the  multipath  affected  signal  and  calculates  the 
off-boresight  angle  to  the  target. 

Secondly,  a  pitch  plane  engagement  model  utilising  the  multipath  affected  returns  from  the  target 
was  developed  and  implemented  in  MATRIXx.  Only  multipath  effects  in  the  vertical  plane  were 
considered.  The  target  is  assumed  to  be  non-manoeuvring  and  travelling  at  constant  speed  and 
height.  It  is  further  assumed  that  the  missile  follows  a  proportional  navigation  guidance  law  and 
travels  at  a  constant  speed. 

The  structure  of  the  report  is  as  follows.  Section  2  discusses  multipath  propagation  in  general,  and 
identifies  factors  that  influence  the  effect,  particularly  in  reference  to  the  reflection  characteristics 
of  the  sea  surface.  Section  3  describes  the  model  used  to  quantify  multipath  propagation  contam¬ 
ination  of  homing  signals  in  a  semi-active  radar  guided  missile.  Section  4  presents  and  discusses 
the  results  obtained  with  the  model.  Finally,  a  brief  summary  of  the  report  conclusions  is  given 
and  the  direction  of  future  investigations  indicated. 


2.  MULTIPATH  PROPAGATION 

The  propagation  of  radar  signals  between  a  source  (the  illuminator),  a  target,  and  a  receiver  (which 
may  or  may  not  be  co-located  with  the  illuminator),  is  influenced  by  several  factors.  These  factors 
include  atmospheric  eflFects  such  as  attenuation  and  multipath  propagation.  The  latter  effect  is  the 
only  one  discussed  here. 

Multipath  is  defined  as  the  propagation  of  electromagnetic  waves  emanating  from  one  point  and 
arriving  at  another  via  different  paths.  In  the  case  of  radar  propagation,  effects  due  to  multipath 
occur  when  two  or  more  electromagnetic  waves  emanating  from  a  single  source,  propagate  along 
different  paths  to  the  same  receiver. 
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Figure  1  The  multipath  engagement  scenario  for  a  generic  missile  against  a  low  altitude  target 

Any  study  of  the  multipath  effects  on  a  semi-active  radar  guidance  system  must  first  recognise  that 
this  involves  a  bistatic  radar  wherein  the  transmitter  is  on  the  launch  platform  (a  ship)  and  the 
receiver  is  in  an  airborne  missile.  This  gives  rise  to  differences  in  the  possible  propagation  paths, 
target  cross  section,  doppler  shifts,  and  glint,  as  compared  to  a  monostatic  radar  wherein  the 
transmitter  and  receiver  are  co-located  [1].  For  example,  consider  the  multipath  scenario  depicted 
in  Figure  1.  When  the  target  is  illuminated  by  the  shipboard  tracking  radar,  the  direct  reflection 
of  this  illumination  from  the  target  to  the  missile  occurs  from  a  positive  aspect  angle  relative  to 
the  illumination  axis  and  thus  has  a  different  effective  radar  cross  section  (RCS)  than  that  seen  by 
the  target  illuminator.  Another  target  reflection  which  propagates  toward  the  sea  surface  occurs 
from  a  negative  aspect  angle  with  respect  to  the  illumination  axis.  This  effective  cross  section  is 
usually  different  from  that  of  the  direct  path  and  can  even  be  larger  for  targets  which  have  more 
radar-reflective  lower  structures.  The  bistatic  differences  in  glint  and  doppler  shifts  arise  from  these 
same  differences  in  the  aspect  angles  for  the  direct  and  reflected  paths  versus  the  illumination  axis 
[1]- 

The  strength  of  the  sea  surface  reflected  signal  depends  on  several  factors.  These  include  the 
efficiency  of  the  reflection  process  (wavelength  and  sea  state  dependent)  and  the  polarisation  of  the 
incident  radiation.  Many  systems  operate  at  X-band,  hence  a  wavelength  of  3  cm  (10  GHz)  was 
adopted.  The  effect  of  these  parameters  is  briefly  indicated  below  (sections  2.1  to  2.3),  as  well  as 
a  discussion  of  the  qualitative  effects  of  multipath  returns  on  missile  guidance  (section  2.4). 
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2.1  Polarisation 

In  order  to  minimise  multipath  effects,  many  semi-active  radar  systems  use  vertical  polarisation. 
This  is  because  horizontally  polarised  radar  at  X-band  is  reflected  far  more  efficiently  from  a 
flat  surface  than  vertically  polarised  radar.  Hence  the  latter  is  likely  to  result  in  significantly 
weaker  multipath  signals. 

2.2  Smooth  sea  reflection  coefficient 

The  efficiency  of  the  reflection  process  is  dependent  on  the  grazing  angle,  the  angle  between 
the  incident  radiation  and  the  reflecting  surface,  as  illustrated  in  Figure  1.  On  reflection  from  the 
sea-surface,  the  negative  aspect  angle  wave  will  suffer  both  attenuation  and  a  phase  change.  For 
vertical  polarisation,  this  attenuation  in  magnitude  and  attendant  phase  change  is  described  by 
the  smooth  sea  reflection  coefficient,  also  known  as  the  Fresnel  reflection  coefficient  [4],  namely: 


€csin^  —  —  cos^  ^ 

CcSinV’  +  \/cc  —  cos^  V’ 


(1) 


where  the  complex  dielectric  constant  €c  appearing  in  the  above  expression  is  approximated  by 


Cc  =  e  -  j60A(t,  (2) 

and  e  is  the  sea  dielectric  constant,  A  is  the  radar  wavelength  (m)  and  a  is  the  sea  conductivity 
(mho/m). 

The  magnitude  of  this  reflection  coefficient  as  a  function  of  grazing  angle  has  been  calculated  for 
typical  sea  water  properties  and  reproduced  as  Figure  2.  The  corresponding  phase  for  similar 
grazing  angles  has  been  depicted  in  Figure  3.  From  Figure  2,  it  is  clear  that,  at  the  so-called 
‘Brewster  angle’,  the  reflection  of  vertically  polarised  radar  at  X-band  is  significantly  reduced 
and  consequently  the  indirect  or  reflected  wave  will  be  attenuated  far  more  than  at  other  angles. 

2.3  Sea  state 

Reflection  from  the  sea  surface  may  be  specular  or  diffuse,  or  a  combination.  Specular  reflection 
occurs  from  a  surface  which  is  flat  and  very  smooth.  If  the  surface  is  irregular,  the  multipath 
is  called  diffuse.  In  the  case  where  the  surface  is  smooth,  but  is  perturbed  by  small  scale 
irregularities,  both  specular  and  diffuse  multipath  are  present,  producing,  in  effect,  coherent 
and  incoherent  components  respectively.  For  the  purpose  of  this  report,  only  the  specular 
multipath  return  is  considered.  This  corresponds  to  ideal  conditions  in  which  the  sea  surface  is 
very  smooth  and  is  generally  designated  as  sea  state  zero.  The  effects  of  the  diffuse  components 
are  discussed  in  [3]. 

In  the  case  being  considered  here,  with  low-altitude  targets  approaching  at  shallow  angles  over 
the  reflecting  surface  of  the  sea,  there  may  not  be  enough  separation  in  angle,  range,  or  floppier 
to  resolve  the  direct  and  reflected  waves  from  the  target.  As  the  multipath  is  specular,  this  wifi 
represent  a  situation  equivalent  to  a  number  of  unresolved  targets,  which  are  related  according 
to  the  geometry  and  electrical  characteristics  of  the  reflecting  surface.  This  will  cause  a  major 
problem  for  the  guidance  of  the  missile  as  the  apparent  target  position  will  vary  with  time  with 
the  variation  in  the  direct  and  reflected  signals. 
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SMOOTH  SEA;  FREQ  =  lOGHz;  VERTICAL  POLARIS AHON 


Pigure  2  The  reflection  coefficient  for  reflection  from  the  sea  surface 


SMOOTH  SEA;  FREQ  =  lOGHz;  VERTICAL  POLARISATION 


Figure  3  The  phase  angle  change  on  reflection  from  the  sea  surface 
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2.4  Effect  of  multipath  returns  on  missile  guidance 

The  signal  received  at  the  missile  is  processed  to  determine  its  direction  of  origin.  In  the  presence 
of  multipath,  this  signal  is  a  combination  of  the  direct  target  skin  return  and  a  return  seemingly 
emanating  from  the  target  image  beneath  the  sea  surface.  With  a  monopulse  receiver,  the 
sum  and  difference  voltages  generated  by  this  composite  signal  are  electronically  manipulated  to 
obtain  the  off-boresight  angle  to  the  target,  that  is,  the  angle  between  the  pointing  direction  and 
the  ‘apparent’  direction  to  the  target.  This  angle,  denoted  here  by  Cm,  is  processed  to  provide 
appropriate  guidance  signals  for  the  missile. 

From  Sherman  [5],  the  indicated  angle  to  the  apparent  target  is  proportional  to  the  ratio  of 
difference  and  sum  voltage  gains  and  is  given  by: 


Cm  —  K  Ji,C 


where  K  refers  to  the  slope  of  the  monopulse  response  curve  versus  off-boresight  angle  and  A,  E 
denote  the  diflFerence  and  sum  voltages,  respectively. 

It  can  be  shown  that  for  the  case  under  consideration  here,  Cm  is  well  represented  by  [5]: 

_  (t  +  p{u  +  gj)  cos<l)  4-  U  ..s 

1  -|-  2  p  cos  <f>  + 

where  Ct  and  Cj  are  the  oflF-boresight  angles  to  the  target  and  image  respectively,  </>  is  the  sum 
of  phase  angles  of  the  direct,  (f>d,  and  reflected,  (l>r,  returns,  that  is: 

+  <^r 

and  p  is  the  ratio  of  the  amplitude  of  the  image  signal  with  respect  to  the  target  signal  (see 
section  3.3). 

The  behaviour  predicted  by  this  equation  is  examined  in  the  following  where  the  open  loop  case 
is  considered,  that  is,  the  seeker  head  remains  fixed  and  only  the  open  loop  indicated  angle 
(relative  to  the  seeker  boresight)  varies  as  a  result  of  changes  in  the  relative  phase,  amplitude 
ratio,  and/or  angular  locations  of  the  target.  In  this  case  it  is  assumed  that 

Ct  =  -U 


and  therefore  equation  (4)  reduces  to 


1  -|-  2  p  cos  (f>  +  p"^ 


Equation  (5)  is  represented  graphically  in  Figure  4  for  a  range  of  indicative  values  of  p  from  0.1 
through  to  5. 

Figure  4  shows  that  the  apparent  line-of-sight  to  the  target  can  be  quite  large,  especially  at 
relative  phase  angles  near  180°,  as  well  as  being  a  periodic  function  of  the  relative  phase  between 
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Relative  phase  angle  (degrees) 


Figure  4 


The  variation  of  the  apparent  off-boresight  angle  as  a  function  of  the  relative 
phase  between  direct  and  indirect  target  returns  as  given  by  equation  (5) 


the  returned  signals.  Furthermore,  for  values  of  p  <  1,  the  apparent  target  always  lies  above 
the  sea  surface  (defined  by  =  0),  and  conversely,  for  p  >  1,  the  apparent  target  lies  below 

the  sea  surface.  For  p  near  1,  the  apparent  target  lies  near  the  sea  surface  except  for  relative 
phases  near  180°.  This  latter  point  indicates  that  if  p  is  varying  about  unity,  then  the  angular 
error,  can  fluctuate  wildly  between  two  extremes  which  extend  beyond  the  angle  between 
the  target  and  image.  However,  because  of  the  relatively  long  seeker  time  constant  and  high 
closing  velocities  associated  with  a  missile  engagement,  this  is  a  transient  situation  which  should 
cause  little  degradation  of  the  missile  performance. 

It  should  be  noted  that  this  was  for  an  open  loop  response.  In  the  closed  loop  case  which  applies 
here  (see  Figure  7),  the  seeker  attempts  to  null  the  angle  error  this  having  the  effect  of 
varying  p,  thereby  making  the  analysis  more  complicated  [5].  However,  the  simpler  open  loop 
response  has  been  used  to  indicate  qualitatively  the  factors  which  will  affect  missile  guidance  in 
a  multipath  situation. 


3.  MODELLING  THE  ENGAGEMENT 

In  this  section,  the  components  of  the  model  are  discussed.  These  components  include  the  engage¬ 
ment  geometry,  the  model  of  the  receiving  antenna,  the  guidance  law,  the  reflection  characteristics 
of  the  sea  surface,  and  the  target  model.  Finally,  the  modelling  package  used  for  these  investigations 
is  discussed  including  a  break  down  of  the  major  sections  of  the  model. 
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Figure  5  The  geometry  for  the  multipath  engagement  scenario 

3.1  Engagement  geometry 

To  simplify  the  discussion  which  follows,  a  flat-earth  geometry  is  assumed.  The  geometry  for  the 
engagement  under  consideration  is  shown  in  Figure  5.  The  grazing  angle  of  the  reflected  wave, 
V’,  will  have  a  major  bearing  on  the  reflection  coefficient  and  phase  angle  change  on  reflection 
from  the  sea  surface  (see  Figures  2  and  3).  This  angle  can  be  determined  from  the  geometry  to 
be  given  by: 


where  h^  and  h„i  are  the  heights  of  the  target  and  missile  respectively,  and  R  is  the  length  of 
the  direct  path  between  the  target  and  missile.  It  can  be  similarly  shown  that  the  length  of  the 
reflected  path,  Ri  -|-  R2)  is  given  by: 

Ri  +  Rj  =  (7) 

sin  xp 

The  phase  of  the  signal  incident  at  the  receiver  is  an  important  quantity  for  reconstructing  the 
direction  of  arrival  of  the  return.  From  the  geometry,  the  phase  of  the  direct  return  from  the 
target,  (j)d,  is  given  by: 

,  27rR 

•Pi  =  — 

and  the  phase  of  the  return  from  the  sea  surface,  <f>r,  will  be: 


(8) 
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(j)j. 


2  TT  (Rj  +  R2) 
A 


+  A<?!)(V0 


(9) 


where  A  is  the  wavelength  of  the  illuminating  radar  and  A<f>{'ip)  is  the  phase  change  on  reflection 
for  incidence  angle  ip. 


3.2  Received  signal 


If  we  consider,  in  general,  the  direct  return  to  have  an  amplitude  Ad  and  phase  <pd,  and  the 
indirect  return  to  have  an  amplitude  Ar  and  phase  then  the  two  signals  can  be  represented 
by  the  time  varying  phasors 

Ad  sin(u;t  +  4>d) 


and 


Ar  sin(a;t  +  4>r) 


Therefore,  the  total  received  signal  will  be: 


At  sin(a;t  +  </>/)  =  Ad  sin(a;t  -f  4>d)  +  sin(a;t  +  <j)r)  (10) 

where  At  and  (j)t  are  the  amplitude  and  phase  of  the  total  received  signal.  Expanding  this 
equation  and  equating  terms  gives: 

At  sm<j)t  =  Ad  sm<l>d  +  Ar  sin<;!!>r 


and 


At  cos  4>t  —  Ad  cos  (f)d  +  Ar  cos  <f)r 


Then  in  terms  of  the  direct  and  indirect  signals,  the  received  signal  will  be  characterised  by  an 
ampUtude  and  phase  given  by: 


A?  = 


and 


{Ad  sin  (/id  +  Ar  sin(f)rf  +  {Ad  cos  (l>d  +  Ar  COS^r)^ 


tan  6i 


_  Ad  sin  (/id  +  Ar 


\  A\ 


(11) 

(12) 


3.3  The  antenna  gain  pattern 

The  receiving  antenna  has  a  gain  which  is  dependent  on  the  angle  to  the  boresight  at  which  the 
incident  radiation  arrives,  the  off-boresight  angle,  e.  For  a  generic  monopulse  antenna  pattern, 
the  gain  function  of  the  voltage  sum  pattern,  G(€),  may  be  well-described  by  the  following 
expression : 

G(c)  =  (13) 

X 

where  Ji(a;)  is  a  Bessel  function  of  the  first  kind  of  order  one.  The  gain  function,  as  represented 
by  this  expression  with  x  =  Ottc  is  illustrated  in  Figure  6. 
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NORMALISED  ANTENNA  GAIN  2-J1  (X)/X 


Figure  6  An  illustration  of  the  model  of  the  antenna  beam  pattern  adopted 


The  gain  of  the  antenna  will  act  to  attenuate  the  signal  strength  entering  the  receiver.  If  G(et) 
and  G(e,)  are  the  antenna  gains  in  the  direction  of  the  target  and  image  respectively,  then  the 
relative  amplitude  of  the  image  signal  with  respect  to  that  of  the  target,  p,  is  given  by: 


p  =  Pip2{i’) 


G(et) 


(14) 


where 

Pi  is  the  ratio  of  radar  cross  sections  of  the  image  and  target  reflecting  surfaces  as  viewed 
by  the  missile  receiver,  and 

P2  is  the  magnitude  of  the  Fresnel  reflection  coefficient  dependent  on  the  grazing  angle  tp  as 
given  by  equation  (1). 

If  we  consider  =  0  as  our  reference  phase,  and  p  to  be  the  ratio  of  signal  amplitudes  of  the 
image  to  target  as  defined  by  equation  (14),  then  equations  (11)  and  (12)  can  be  rewritten  as 

=  1  +  +  2pcos0  (15) 


where  (p  =  +  <f>rt  and 


tan  (pt 


p  sincp 
1  +  p  cos<p 


(16) 
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Equations  (15)  and  (16)  can  be  used  to  calculate  the  phase  and  amplitude  of  the  total  received 
signal  from  the  target. 


3.4  Generic  missile  guidance 

A  basic  pitch  plane  trajectory  model  is  used  for  the  missile.  A  brief  description  of  the  model 
follows. 

The  inertial  angles  relating  the  angular  positions  of  the  target,  the  image,  0i,  and  the  pointing 
direction  of  the  seeker  head,  Bg,  are  used  to  determine  the  off-boresight  angles  to  the  target  and 
image,  ct  and  e,-.  These  are  used  with  the  received  signal  represented  by  p  and  <j),  and  processed 
in  accordance  with  equation  (4)  to  determine  the  apparent  boresight  error  Cm-  This  is  then  used 
for  missile  guidance  based  on  a  proportional  navigation  (PN)  guidance  algorithm.  In  this  case, 
the  turning  rate  of  the  missile,  7,  is  proportional  to  the  turning  rate  of  the  missile-to-target 
line-of-sight  (LOS),  &,  that  is: 

7  =  ka  (17) 

where  k  is  the  proportionality  constant  called  the  navigation  ratio,  or  PN  gain. 

The  line-of-sight  rate  is  generated  by  the  seeker  head  tracking  the  target  and  attempting  to  keep 
the  target  on  boresight.  Any  off-boresight  deviation,  Cmi  results  in  a  demand  to  move  the  seeker 
head  to  null  the  deviation.  The  off-boresight  angle  is  used  to  reconstruct  the  line-of-sight  rate, 
(T,  using  the  algorithm  illustrated  by  Figure  7,  where  the  seeker  boresight  error  is  combined  with 
the  integrated  rate  gyro  signal  (see  reference  [6]).  The  response  of  the  seeker  gyro  is  modelled 
as  a  first  order  lag  with  a  time  constant  T*. 

The  line-of-sight  rate  is  then  used  to  demand  the  called-for- acceleration,  Nc  (or  A^em  iii  Fig¬ 
ure  7),  given  by: 

Ne  =  N'&y,  (18) 

where  Vg  is  the  relative  closing  velocity  between  the  missile  and  the  target,  and  N'  is  the  effective 
guidance  gain  with  a  value  in  the  range  2-5.  A  value  of  4  is  adopted  here.  This  demand,  with 
g-bias,  is  fed  to  the  autopilot  which  is  assumed  to  respond  as  a  first  order  lag  with  a  time 
constant  of  To,  resulting  in  the  achieved  acceleration. 

3.5  Target  model 

The  target  model  adopted  in  these  preliminary  investigations  is  of  a  non-manoeuvring  target 
travelling  at  constant  speed  and  at  a  constant  altitude.  The  ratio  of  the  target  radar  cross  section 
(RCS)  between  the  positive  and  negative  aspect  angles  is  assumed  to  be  a  variable  input.  In  this 
context,  the  positive  aspect  angle  is  defined  as  the  angle  from  the  illuminator  beam  to  the  target 
echo  towards  the  missile  measured  in  a  clockwise  direction  (see  Figure  1).  Likewise,  the  negative 
aspect  angle  is  defined  as  the  angle  from  the  illuminator  beam  to  the  target  echo  towards  the 
sea  surfcLce  and  measured  in  a  counterclockwise  direction. 

Consequently,  the  RCS  Ratio  is  defined,  as  given  in  reference  [7],  in  terms  of  the  target’s  average 
radar  cross  section  in  the  direction  of  the  line-of-sight  of  the  missile,  cta,  and  its  average  radar 
cross  section  in  the  direction  of  the  specular  reflection  point  on  the  sea  surface,  ctb,  that  is: 


RCS  Ratio 


(19) 


In  addition,  as  mentioned  previously,  the  simulation  program  has  the  facility  to  accommodate 
the  target’s  bistatic  RCS  profile,  should  such  information  be  known  for  a  specific  target. 
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Figure  7  The  guidance  algorithm  modelled  including  line-of-sight  reconstruction 

3.6  Sea  surface  reflection  characteristics 

The  properties  associated  with  the  reflection  of  radar  waves  by  the  sea  surface  are  dependent  on 
the  angle  of  incidence  of  the  radiation,  its  polarisation,  and  the  sea  state.  In  the  case  considered 
here  of  vertical  polarisation  and  a  flat  sea  surface,  the  attenuation  and  phase  change  on  reflection 
as  a  function  of  the  grazing  angle  are  as  given  in  Figures  2  and  3. 

3.7  Modelling  package 

The  model  is  developed  using  the  MATRIXx  simulation  package.  MATRIXx  is  a  powerful, 
programmable,  matrix  calculator  which  offers  a  graphics  capability.  Within  MATRIXx  the  user 
is  able  to  access  the  SystemBuild  environment  where  block  diagrams  are  generated  and  the 
simulation  environment  provided. 

The  MATRIXx  model  to  represent  the  effects  of  multipath  interference  is  represented  in  Figure  8 
by  its  SystemBuild  block  diagram.  This  shows  the  different  blocks  which  make  up  the  flow  of  the 
model.  There  are  six  main  ‘Superblocks’  which  essentially  equate  to  functions  or  subroutines  in 
other  programming  languages.  There  are  also  some  superblocks  within  these  superblocks.  The 
complete  listing  is  not  given,  but  the  calculations  performed  by  each  superblock  are  described 
in  Table  1.  The  multipath  model  is  part  of  the  overall  engagement  model  and  uses  the  position 
of  the  missile  (XM,ZM)  and  the  position  of  the  target  (XT,ZT)  to  determine  the  ‘apparent 
position’  of  the  target  in  the  presence  of  multipath  returns.  This  apparent  position  is  used  by 
the  missile  model  (Figure  7)  to  generate  guidance  commands.  This  will  allow  an  assessment  of 
the  response  of  the  missile  to  target  position  fluctuations  which  occur  in  an  environment  where 
multiple  returns  from  a  target  are  possible.  This  is  examined  in  the  following  section. 
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Table  1  List  of  superblocks  in  program  structure 


Superblock 

Function 

Flat  Earth  Geometry 

In  this  superblock,  the  geometrical  parameters  of  the  engage¬ 
ment  are  calculated.  This  includes  the  range  of  the  missile  to 
the  target  and  image,  and  the  lookdown  angles  to  the  target 
and  image  relative  to  the  missile  centerline. 

— 

Reflection  Coefficient 

This  superblock  determines  the  coefficients  for  amplitude  atten¬ 
uation  and  phase  change  of  the  target  signal  reflected  from  the 
sea  surface  along  the  line-of-sight  to  the  missile  receiver. 

Monopulse  Processor 

The  superblock  determines  the  off-boresight  angles  to  the  target 
(eps  t)  and  the  target  image  (eps  i)  from  the  geometry,  and  the 
off-boresight  angle  to  the  apparent  target  (eps  m)  using  phase 
monopulse  techniques  to  process  the  multipath  signals  received 
by  the  missile  from  the  ‘target’. 

LOS  Rate  Reconstruction 

Using  the  boresight  error  to  the  ‘apparent  target’,  the  line-of- 
sight  rate  is  calculated  assuming  a  seeker  head  time  constant  of 
0.2  seconds. 

Seeker  Head  Angle 

In  this  superblock,  the  seeker  head  angle  is  determined  from  the 
seeker  head  rate  and  line-of-sight  rate. 

Apparent  Target  Height 

Calculates  the  apparent  target  position  from  the  signal  received 
by  the  missile  receiver.  This  uses  the  apparent  target  position, 
the  seeker  head  angle  and  the  position  of  the  missile  relative  to 
the  true  target  position. 

4.  RESULTS  AND  DISCUSSION 

The  engagement  model  was  run  for  a  scenario  chosen  as  representative  of  the  final  stages  of  a  typical 
low-altitude  target  engagement.  The  important  quantities  to  be  extracted  from  the  simulation 
include  those  factors  that  are  indicative  of  missile  performance.  For  the  multipath  environment 
case,  a  significant  factor  affecting  missile  performance  is  the  apparent  target  position  -  the  target 
position  as  seen  by  the  seeker  -  since  missile  guidance  ultimately  depends  on  this. 

The  nominal  parameters  of  the  engagement  scenario  considered  in  this  analysis  are  given  in  Table 
2.  It  should  be  noted  here  that  the  scenario  is  assumed  to  begin  in  the  terminal  phase  of  the 
engagement,  therefore  initial  parameters  relate  to  the  beginning  of  this  phase,  not  missile  launch 
conditions. 

For  the  initial  conditions  given  in  Table  2,  several  cases  of  target  height  and  sea  state  have  been 
considered.  Target  heights  of  5m,  15m  and  30m  have  been  chosen  for  simulation  purposes  as  well 
as  sea  states  of  zero  and  three.  The  effect  of  sea  state  is  represented  by  replacing  p  in  equation  (14) 
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Reflection  Coefficient 


Figure  8  SystemBuild  block  diagram  of  the  multipath  model 
Table  2  Nominal  parameters  for  the  simulated  engagement 


Parameter 

Initial  Value 

Missile  height 

Missile  to  target  range 
Missile  velocity 

Target  velocity 

Wavelength 

Dive  Angle 

Seeker  lookdown  angle,  Og 
RMS  noise 

RCS  ratio  {(TaI(Tb) 

Time  constants  T^ 

Time  constant  T), 

585m 

5400m 

700  m/s 

300  m/s 
0.03m 

9° 

6'’ 

0.0“ 

0  dB 

0.2  sec 

0.25  sec 

by  Pss,  where 


Pss  —  P  Pc 


and  Pc  is  the  coherent  scattering  co-efficient  and  has  a  value  which  varies 
for  sea  state  zero  and  less  than  0.1  for  sea  state  three  [4]. 


non-linearly,  being 


(20) 

1 

^  i 


The  effect  of  multipath  returns  from  the  target  on  the  apparent  target  height,  boresight  error 
variations  and  the  effect  of  sea  state  are  presented  in  figures  9  to  12  and  discussed  below.  It  should 
be  noted  here  that  Cm  and  the  apparent  target  height  are  intimately  related. 
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Figure  9  Apparent  target  height  variations  for  true  heights  of  5m,  15m  and  30m  for  sea  state  zero 


4.1  Apparent  Target  Height 


Figure  9  shows  the  apparent  target  height  variations  for  true  target  heights  of  5m,  15m  and  30m 
and  sea  state  zero.  The  time  axis  indicates  the  time  from  the  commencement  of  the  terminal 
phase  and,  given  a  closing  velocity  of  about  1000  m/s,  the  range  between  the  missile  and  target 
can  be  deduced  using  the  intercept  time  of  approximately  5.45  s.  The  deviation  of  the  apparent 
target  height  from  the  true  target  height  is  due  to  multipath  returns.  On  average,  the  behaviour 
indicates  that  the  apparent  height  closely  follows  the  true  height,  with  multipath  returns  having 
their  greatest  effect  in  the  final  2  seconds  before  intercept.  Variations  of  the  apparent  height 
from  the  true  height  can  be  up  to  about  30%  in  each  case.  The  frequency  of  the  apparent 
target  height  fluctuations  reflect  the  rate  of  change  of  the  path  length  difference  between  the 
target  and  target  image  returns.  In  the  final  stages  of  the  engagement,  where  the  angular 
displacement  between  the  target  and  image  increases  more  rapidly  (see  Figure  10),  the  apparent 
height  variation  is  seen  to  be  the  greatest.  However,  at  no  stage  during  the  simulations  did  the 
apparent  height  converge  towards  the  image  height,  namely,  -5m,  -15m  and  -30m,  respectively. 
This  stems  directly  from  the  magnitude  of  the  total  reflection  coefficient,  which,  in  this  case  is 
sa  0.2,  and  therefore  significantly  reduces  the  strength  of  the  target  image  return. 


4.2  Boresight  error  variations 

Figure  10  presents  the  off-boresight  angles,  or  boresight  errors  (BSE),  to  the  target,  the  target 
image,  Ci,  and  the  apparent  target,  Cm-  For  each  of  the  cases  considered  here,  the  apparent  target 
boresight  error  effectively  follows  the  target  boresight  error.  As  for  the  apparent  target  height, 
the  deviation  of  from  €t  is  due  to  multipath  returns,  and  the  frequency  of  the  variation  is 
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Figure  10  The  apparent  target  boresight  error  ( - ),  the  true  target  boresight  error  ( - 

— ),  and  the  target  image  boresight  error  ( - )  for  target  heights  of  5m,  15m  and  30m  and 

sea  state  zero 
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again  due  to  the  rate  of  change  of  the  path  length  difference  between  the  target  and  target 
image  returns.  Again,  in  the  final  stages  of  the  engagement  the  variation  was  greatest,  but  on 
average,  was  still  centred  on  the  true  target  position.  As  in  the  apparent  height  consideration, 
at  no  stage  did  the  apparent  target  position  tend  to  the  image  position. 

4.3  Sea  state 

In  considering  the  effect  of  sea  state,  entering  through  equation  (20),  a  target  height  of  15m 
was  considered,  and  the  apparent  target  height  variation  and  the  relative  power  of  the  returned 
signal  were  examined  for  sea  states  zero  and  three  as  shown  in  Figures  11  and  12.  A  relative 
power  of  0  dB  corresponds  to  the  direct  target  return  and  variations  from  this  are  due  to  the 
interference  effects  of  the  multipath  returns. 

The  effect  of  sea  state  is  clearly  evident.  A  higher  sea  state  reduces  the  coherent  (specular) 
multipath  returns  and  this  is  reflected  in  Figures  11  and  12.  The  coherent  scattering  coefficient, 
pc  in  equation  (20),  has  a  value  of  w  1  for  sea  state  zero  and  reduces  significantly  to  less  than 
0.1  for  sea  state  three.  Hence  the  observed  effects. 


5.  CONCLUSIONS 

In  this  paper  we  have  presented  a  model  for  assessing  the  effect  of  multipath  returns  on  semi-active 
missile  guidance.  The  model  may  be  used  to  assess  the  performance  of  particular  semi-active  radar 
homing  missiles  when  operating  over  a  sea  surface.  At  present  we  have  only  considered  coherent 
(specular)  multipath,  but  the  model  is  being  extended  to  include  diffuse  multipath  returns  [3]. 

The  results  obtained  above  indicate  that  for  the  initial  conditions  considered  as  representative  of  a 
low-altitude  engagement,  an  examination  of  apparent  target  height  (and  effectively  apparent  target 
boresight  angle)  indicates  that  although  multipath  returns  are  impacting  on  the  predicted  target 
position,  average  behaviour  tends  to  the  true  target  position.  Any  significant  deviations  from  the 
true  target  position  are  extremely  transient  and  are  damped  by  the  missile  response  time.  Therefore 
they  do  not  cause  any  adverse  guidance  problems  for  the  missile. 

A  detailed  assessment  of  the  final  intercept  geometry  (in  terms  of  miss  distance)  is  not  justified  with 
the  simple  missile  model  adopted  here,  as  only  qualitative  effects  were  sought.  A  more  quantitative 
analysis  of  the  intercept  geometry  has  been  pursued  using  a  six  degree-of-freedom  missile  model 
with  the  incorporation  of  diffuse  multipath  returns  (see  reference  [8]  for  a  more  detailed  discussion). 
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